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CO ■ Abstract 

o 

""^ ' The behaviour of a massless Dirac field on a general spacetime background representing 

two colliding gravitational plane waves is discussed in the Newman-Penrose formalism. 



The geometrical properties of the neutrino current are analysed and explicit results are 



given for the special Ferrari-Ibahez solution. 



^ ■ 1 Introduction 

h ; 

Exact solutions of the Einstein equations representing colliding gravitational plane waves have 
been discussed extensively in the literature and the status of research on this topic is still 
best summarized in a monograph by Griffith 1 published in 1991. 

The interest in such solutions is mainly due to the possibility they offer to understand 
better some nonlinear features of the gravitational interaction. In general, the spacetime 
geometry associated with two colliding gravitational plane waves is very rich: for example 
either a spacetime singularity [2] or a Killing-Cauchy horizonjSl EJ O2J can result from the 
nonlinear wave interaction. In general, such spacetimes contain four regions: a Minkowki 
region, representing the initially flat situation before the passage of the two oppositely directed 
plane waves, two Petrov type N regions, corresponding to the waves before the interaction, 
and an interaction region, generally of Petrov type I. Two commuting spacelike Killing vectors 
are always present, associated with the plane symmetry assumed for the two colliding waves. 

Here we consider a massless Dirac field (neutrino) interacting with the colliding plane 
waves. Explicit results are given for the degenerate solutions found by Ferrari and Ibanez^], in 
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which the interaction region of the colliding waves propagating along the common z direction 
is of Petrov type D, with a Killing- Cauchy horizon formed at a finite distance from the collision 
plane, after a finite time from the instant of collision. For this spacetime a careful analysis of 
timelike geodesies has only recently been performed 6 . 

Some special features of the neutrino current are discussed in detail, and previous work 
of Dorca and Verdaguer[71 |S] and YurtseverjH] valid for a scalar field interacting with gravi- 
tational waves is extended here to spin 1/2 massless particles. 



2 Colliding waves: the Dirac equation 

The most general form of the spacetime metric representing two colliding waves with parallel 
polarization can be written as 

ds 2 = 2g 12 dxi dx 2 - 533 (dx 3 ) 2 - g u (dx 4 ) 2 , (1) 

where g\ 2 , 533, 544 are real functions of the null coordinates (xi, x 2 ). The complete spacetime 
description of the collision must patch together two single-wave regions (the approaching 
waves) plus a portion of flat spacetime corresponding to the situation before the passage of 
the waves and the nontrivial interaction region where the metric depends on both x\ and 
x 2 . We distinguish the single- wave metrics by their dependence on a single null coordinate 
(,9ij = 9iji x l) f° r t ne progressive wave and g^ = gij(x 2 ) for the regressive one), while in flat 
spacetime these components are constants, g\ 2 = 2, 333 = 544 = 1. Khan and Penrose [2] 
identified a standard procedure to extend continuously the metric over all the various regions 
by a proper use of Heaviside step functions. We will discuss the details of this extended 
spacetime in the case of the explicit solution of Ferrari-Ibahez in the next section. 

A wide variety of exact solutions can be found in the literature, the analysis of which is 
extremely simplified using the Newman-Penrose formalism ^U]. Here the following NP null 
frame can be introduced 



I = d X2 , n = I/512 d Xl , m = —= 



1 „ 1 



9x3 ^ 1 



/933 yj9U 

It is convenient to introduce the following two quantities 



(2) 



S = VWW , $ = vW(2S44) (3) 
which play a role in a conformally rescaled form of the metric Q 

ds 2 = ^ [4ctei dx 2 - H~ 2 (dx 3 ) 2 - <5>~ 2 (dx A ) 2 ] . (4) 

The Dirac equation for massless spin 1/2 particles then reduces to 

(D + e - p) F x + (5* + 7T - a) F 2 = , 

(S + p - t) F 1 + (A + n - 7) F 2 = , (5) 

where the notation and conventions for the NP formalism follow those of Chandr asekhar [TT] . 
A much simplified form of (JSJ) results from rescaling the spin wave functions Fx and F 2 as 
follows 

Fx - G?33 94 4 r 1/4 /v^#i = ("5)~ 1/4 H, , 

F 2 -» (<733 544)- 1/4 ^2 • (6) 
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According to the Khan-Penrose procedure 2, we must first examine the Dirac equation for 
each of the four regions. 



2.1 Case 1: g, Lj = g^Xi, x 2 ) 

The NP spin coefficients have a very simple and symmetric form in the general case 



P=—g8x2 lo g (533 544) 

\ = -^-d Xl log 

4^12 944 

e = 2 5*2 (log 012) 
while the Weyl scalars are 

^ = Da- 2a(p + e) = 



O- = "J ^2 log( , 

4 544 



1 



9xi lo g (#33 544) 



4^12 

a = /3 = 7 = K = i/ = 7r = T = 



[5^2 !og (512 533)] 2 - [9x2 lo g (#12 S'44)] 5 



933 9l 2 44 



533 



544 



^1 = -5e 
Dfj,- Ae + 2(/ze- Act) 



1 



(2 d xi log 012 5z 2 log 512 



12512 

~2 9xi log ^ 533 944 ^ 5x2 log ( Q33 9AA ^ 

_2 9x 2 512 5 Xl <9x 2 033 5a:i 5^2 544 



512 



533 



544 



$3 

^4 = - A A - 2\fi 



[5xi lo g (512 533)] 2 - [5xi lo g (512 044)]' 



5xi #33 9 X1 044 



533 



544 



(7) 



(8) 



Since the metric is independent of (2:3, £4) one can seek solutions (normal modes) of the 
form 

Hi = e i( K ^3+K 4XA ) ^ ^ = e i(K 3 x 3 +K 4 x 4 ) ^ ^ ^ 

where K 3 and K4 are real constants and hi, h 2 will depend on them. The Dirac equation 
then becomes 



dhi{xi, x 2 ) 

9x 2 
9h 2 (x 1 , x 2 ) 
dxi 



(K 4 ^-iK 3 E)h 2 (x 1 , x 2 ) , 
(-K A $ -iK 3 E)h 1 (x 1 , x 2 ) , 



and the general solution can be given as superposition of modes 



Hi = / dK 3 dK 4 Hi , H 2 = dK 3 dK 4 H 2 . 



(10) 



(11) 
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2.2 Case 2: g tj = ftj(x x ) 

In this case, only A and \i remain nonzero among the spin coefficients ([7)l. Since the metric 
is independent of (x2, £3, 24) one can seek solutions of the form 

Hl = e i(K 2 x 2 +K 3 x 3 +K 4 x 4 ) h ^ Xi) ^ ^ = e i(X 2 x 2 +^ 3 X3+X 4 z 4 ) ^ ( 12 ) 

where X2, -K3 and K 4 are real constants and h\ and /12 will also depend on them. The Dirac 
equation then becomes 

hi{xx) = — (-K 3 3-iK A $)h 2 (x 1 ) , 

h' 2 (xi) = ^-[(K 3 E) 2 + (K 4 ^) 2 }h 2 (xi) , (13) 
and the general solution can be given as superposition of these modes 

U x = J dK 2 dK 3 dK 4 Hi , H 2 = j dK 2 dK 3 dK 4 H 2 . (14) 

2.3 Case 3: = gij{x 2 ) 

Here the non-zero spin coefficients are p, a and e. Since the metric is independent of 
(xi, x 3 , X4) one can seek solutions of the form 

Hl = e i(K 1Xl +K 3 x 3 +K 4 x 4 ) ^ ^ = e i(X 1 x 1 +iC 3 X3+X 4 X4) ? ( 15 ) 

where K\, K 3 and K 4 are real constants and as in the previous cases hi and h 2 will also 
depend on them. The Dirac equation then becomes: 

h[(x 2 ) = ^-[(K 3 E) 2 + (K 4 ^>) 2 ]hi(x 2 ) , 

h 2 (x 2 ) = ±- (-K 3 E + iK 4 ^)h 1 {x 2 ) , (16) 
and the general solution can be given as superposition of these modes: 

Hi = J dKi dK 3 dK 4 Hi , H 2 = j dK x dK 3 dK 4 H 2 . (17) 

2.4 Case 4: g^ = (positive) constants 

For a metric with constant components the spin coefficients are obviously all zero. We assume 
without any loss of generality: 512 = 2 and 533 = 544 = 1. Since the metric is constant one 
can seek solutions of the form 

_ e 'i (Ki X1+K2 X2+K3 X3+K4 x 4 ) ^ jj^ _ g i 11+^2 12+^3 13+^4 14) ^ 2 /^g\ 

with hi and h 2 also constants. The Dirac equation reduces to 

K 2 hi + (K 3 + i K 4 ) h 2 = , 

(K 3 -iK i )h 1 +K 1 h 2 = 0, (19) 
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which implies (by setting to zero the determinant of the coefficient matrix of the linear system 

K x K 2 = Ki + Kl . (20) 

In the next section we will specialize these cases to the horizon-forming Ferrari-Ibahez 
metric and discuss their solutions. 



3 The horizon-forming Ferrari-Ibahez metric 

In 1987 Ferrari and IbahezjUEj found a type D solution of the Einstein equations that can 
be interpreted as describing the collision of two linearly polarized gravitational plane waves 
propagating along a common direction z in opposite senses and developing a non-singular 
Killing- Cauchy horizon when they collide. Using the standard coordinates (t, z, x, y), with 
an appropriate choice of the amplitude parameters associated with the strength of the waves, 
the metric takes the form: 

ds 2 = (1 + sin t) 2 {dt 2 -dz 2 )- 1 ~ smt dx 2 _ cos 2 z (i + S i n t )2 dy 2 _ ( 2 1) 

1 + sin t 

The interaction region where this form of the metric is valid (designated as "Region I" , as 
in most of the literature [7]) is represented in the (t, z) diagram by a triangle whose vertex 
(representing the initial event of collision) can be identified with the origin of the coordinate 
system; the horizon is mapped onto the base of the shaded triangle in Fig. 1. 

In order to describe the larger spacetime of which this is only one region, one must 
introduce the two null coordinates 

u = (t - z)/2, v = (t + z)/2 <{=>- t = u + v, z = v-u , (22) 

in terms of which the metric ()21|) takes the form 

ds 2 = 4 [1 + sin(-u + v)] 2 du dv 

dx 2 - cos 2 (ti - v)[l + sin(u + v)] 2 dy 2 . (23) 



1 + sin(n + v) 

Following Khan- Penrose % one can easily extend the formula for the metric from the 
interaction region to the remaining parts of the spacetime representing the single wave zones 
and the flat spacetime zone before the waves arrive. The interaction region corresponds to 
the triangular region in the (u, v) plane bounded by the lines u = 0, v = and u + v = ir/2. 
One need only make the following substitutions in ((2 



u — ► u H{u) v — > v H(v) 

which give rise to the four regions 

u>0, v>0,u + v<tt/2 Region I Interaction region 

< u < tt/2 , v < Region II Single u-wave region 

n<0, < v < it/2 Region III Single v-wave region 

u < , v < Region IV Flat space 



(24) 



shown in fig. 1. In this way the extended metric in general is C° (but not C 1 ) along the null 
boundaries u = and v = 0. It is worth noting that certain calculations are more easily done 
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Figure 1: The null coordinates and the different regions they induce. 



in one or the other of these two sets of coordinates, so we will switch back and forth between 
them as needed. 

We can apply the results of the previous section using the coordinates (u, v, x, y) and the 
extended metric ((23J - (JUJ). 

3.1 Region I: wave interaction 

The metric 

9 i2 = 2 [1 + sin(u + v)} 2 , 
1 — sin(-u + v) 
1 + sm(u + v) 

944 = cos 2 (f — u) [1 + sin(u + v)] 2 . (25) 

is of Petrov type D, the NP frame © is not principal (i.e. not aligned along the two repeated 
principal null directions of the spacetime) and the Dirac equation (|10[) becomes 



dh\(u, v) ( K y [1 + sin(it + v)] 



iK x — ■ — r h 2 (u, v) , 



dv \cos(v — u) cos(n + v 

dh,2(u, v) ( K y [1 + sin(n + v)] 



iK x v - \ / J hi(u, v) , (26) 



du \ cos(v — u) cos(n + v) 

where a convenient notational change has been made with respect to the previous section, 
denoting #1,2,3,4 as K UtV)X) y. 

To exploit the symmetries of (|26p we again transform the Dirac variables 

hi(u, v) = si(u, v) + s 2 (u, v) , 

h 2 (u, v) = si(u, v) - s 2 {u, v) , (27) 
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and then we switch back to the spacetime coordinates (t, z, x, y), introducing the abbreviated 
notation 



a(t) 



(1 + sinty 
cost 



f3(z) 



1 



COS z 



Equations (|26|) then become 



dsi(t, z) ds 2 {t, z) 



dt 



dz 



dsi(t, z) ds 2 (t, z) 



dz 



dt 



-i K x a(t) si (t, z) - K y j3(z) s 2 (t, z) , 
K y P(z)si(t, z) + iK x a(t)s 2 (t, z) . 



(28) 



These equations can be separated by assuming s\(t, z) = Ti(t)Z\(z) and s 2 (t, z) = T 2 (t)Z 2 (z). 
In fact dividing the first by T 2 (t) Z\(z) and the second by T\(t) Z 2 (z), we obtain 



T[{t) , Z' 2 {z 



+ 



+ iK x a{t) 



+ 



Ti(t) Z 2 (z) 



i K x a(t) 



Ti(t) 

r a (t) 

Ti(t) 



K y (3(z) 



Z2(Z) 

Zi(z) 
Zi(z) 

Z2(Z) 



(29) 



and hence 



T[(t) + iif x a(t)Ti(t) 
T^(t)-iK x a(t)T 2 (t) 
Z[(z) - K y P(z) Z x {z) 
Z' 2 {z) + K y P(z) Z 2 (z) 



-KT 2 {t) , 
-JTx(t) , 
JZ 2 (z) , 
KZ^z) , 



(30) 



where J and K are arbitrary constants. By differentiating these equations the functions T, Z 
are seen to satisfy a system of four completely separated second-order differential equations 
of the following form: 



df{x) 
dx 2 



[a 2 A{xf +aA'(x) + JK]f(x) = , 



(31) 



with the respective sets of values for f(x), a, A: 



( Tt 

T 2 
Zi 

{ z 2 



iK x 
iK x 
K y 



A 



f «(*) 
a(t) 

0{z) 

I 0(z) 



(32) 



With the restriction {J K = — (§) 2 , n G N}, the (|31jl for Z±, Z 2 can be integrated analyti- 
cally. In fact with the variable change z = ir/2 — 9 and the rescaling Z\ 2 (6) = Vsin# @i, 2 (0) 
the resulting equation for ®i t2 (9) becomes 



e"(6) + cot6Q'{9) + 



2 s m cos 9 — m? — s 2 



sin 



+ 1(1 + 1) 



Q(9) = 



(33) 



where J K = —[1(1 + 1) + s 2 ] and s = +1/2 must hold in order to have regular solutions, with 
Q\(0) corresponding to s = —1/2, m = K y (or s = 1/2, m = —K y ) and Q 2 (9) corresponding 
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to s = 1/2, m = K y (or s = — 1/2, m = —K y ). A class of solutions of (|33|) is represented by the 
spin-weighted spherical harmonics S Y, v (6 1 0) (first introduced by Newman and Penrose [10] 
as a natural extension of the scalar harmonics to the Hilbert space of the square-integrable 
spin-s functions defined on the 2-dimensional unit sphere). An explicit definition of these 
functions, together with a list of their main properties and features is given in the Appendix. 
The corresponding solutions of Eqs. (|3flH32]) can then be denoted by 

T x (t) = Tf*(u + v), 

T 2 (t) = T j- K *(u + v) 

Z\(z) = ^/cos(v — u) -i/ 2 Y l Kv (n/2 + u — V, 0), 

Z 2 (z) = Vco^ - u) xj2 Y* v (vr/2 + u-v,0) (34) 

where indices 1, 2 correspond to the signs ± respectively. The general solution will then be a 
superposition of such modes 

oo 

h t (u,v) = yfiM^uj C l [Tf*{u + v) 1/2 Y- Ky (ir/2 + u-v,0) + 

l = \Ky\ 

TJ Kx (u + v) 1/2 Y t Ky (vr/2 + u-v, 0)] (35) 

with an analogous expression for h 2 , where but the functions T^ Kx {u + v) can only be given 
numerically. 



3.1.1 Teukolsky Master Equation 



We notice here that as this portion of spacetime is Petrov type D actually one can write, 
following notation and conventions of Teukolski|12j. a Master Equation for any spin s = 
0, ±1/2, ±1, ±3/2, ±2 perturbations of the form 



if>(t, z, x, y) = e iK * x e iR y y T(t) Z{z) 



(36) 



of this background. To do this it is convenient to adopt a (Kinnersley) principal null tetrad, 
different from the one introduced above for the Dirac field, aligned with the two repeated 
principal null directions of the spacetime 



I 



n 



m 



1 



2 cost 
cos t 



(1 + sin£) 2 
1 

\/2(l + sint) 



cos t 
dt - d x , 



d z + 



-d„ 



cos z 



so that one gets two decoupled equations, one for Z(z) 

d 2 Z{z) dZ{z) 



where 



dz 2 



VJz) 



tan z- 



K 2 



dz 



V z (z)Z(z) = , 



2K v s sin z + s z 



cos 2 z 



(37) 

(38) 
(39) 
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with £1 a separation constant, and one for T(t) 



^ - (1 + 2s) tan - V t (t)T(t) = , (40) 

where 

V t (t) = n + s(s + 1) - Kl (1 + ^ - 2 iK x s ( sint + 2^±ff . (41) 

COS z t \ COS z I / 

Again, by introducing the new variable 9 = tt/2 + z Eq. (jHEJ takes the same form of Eq. (j33j) . 
and with = — Z (Z — |— 1 ) , with l,m = K y integers or half-integers and s = 0, ±1/2, ±1, ±3/2, ±2, 
its solutions are the spin-weighted spherical harmonics, with spin-weight s, generalizing the 
results already given for the Dirac field. 

To accomplish the task of extending this description all over the spacetime it is obviously 
necessary to match the equations for the various spin fields also in the other two regions of 
type N and in the Minkowski portion. This can be done for spin s = 1 in complete analogy to 
what we have done for the s = 1/2 case. The gravitational case s = 2 is to be more carefully 
handled because of the nonvanishing radiative contribution of the background gravitational 
field in the two N regions. However this analysis goes beyond the aims of the present paper 
and it will be studied subsequently. 

3.2 Region II: single «-wave 

The metric 

512 = 2 (1 + sinn) 2 , 
1 — sin u 

533 : 



1 + sin u ' 

544 = cos 2 (n) (1 + sinu) 2 , (42) 

is of Petrov type N, the NP frame © is principal and the Dirac equation ()13j) gives 

-K x (1 + sinu) 2 - iK y 

hi(u) = y - h 2 {u) , 

K v cosn 



= i ~ K 2 L h 2 (u). (43) 

K„ cos z u 



The ordinary differential equation for fi2{u) can be solved exactly 

h 2 (u) = -V2N ^ ek l K * A ^ +K2 y ta H , (44) 

K x + i K y 

where TV is a normalization constant and the function A(x) is defined by 

.. . (sinx + 8)(cos 2 x + 2) 15 in fAr . 

A(x) = — + 7 tan x x - 12 . (45) 

v 1 2 cos x 2 v ' 
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3.3 Region III: single n-wave 

The metric 



512 = 2(l+sinv) 2 , 

1 — sin v 

533 = ^— — : j 

1 + sm v 

544 = cos 2 (v) (1 + sin-u) 2 (46) 
is again of Petrov type N, the NP frame © is principal and the Dirac equation (jlfij) is: 

h^v) = i hi{v) , 

K u cos z v 

h 2 (v) = y - h^v) . (47) 

K u cosv 

The solution for h\ (v) is 

ht (v) = V2 N e^t [*% A ^+ K v tan »] , (48) 

where N is a normalization constant and A(y) is defined in (|45() . Here one has the same 
problems matching the solution to region I as for region II. 

3.4 Region IV: flat spacetime 

The metric 

512 = 2, (/33 = 544 = 1 . (49) 

is the Minkowski metric. The relation between the K-constants that comes from the Dirac 
equation (f20|) is 

K u K v = K 2 X + Ky (50) 

(which corresponds to energy-momentum conservation for the test particle^)]). 
A straightforward solution for the amplitudes of the wave functions is: 

hl = V2N, h2 = -^ N -J±- = -^ N >S^. (51 ) 

This list of the various cases concludes the analysis of the Dirac equation. We have seen 
that an explicit solution can be given completely in the two type N regions and of course 
in the flat Minkowski region. In the type D region the Dirac equation can be separated in 
its (t , z) dependence, but only the separated z equation can be solved exactly in terms of 
spin- weighted spherical harmonics. 

3.5 Tracking the Dirac field 

It is now clear how to obtain the formal solution of the Dirac equation in the four spacetime 
regions. However, since the general solution in region I can only be given numerically, the 
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problem of (continuously) matching at the boundaries I — II and / — III can only be ap- 
proached numerically. To be more explicit, consider the example of a transition from region 
I to region II for h\. Approaching the it axis (v — » 0) one has 

oo 

lim hi (u,v) = v/cosn Ci[Tf"(u) 1/2 Yf Kv (n/2 + u, 0) + 

l = \Ky\ 

T^{u) l/2 Y^^/2 + u,Q)] , (52) 

and this limit must coincide with the solution hi(u) of Eq. (|43|) in region II, so that it can 
be seen as a series representation for h\{u) in terms of spin- weighted spherical harmonics. 
But the coefficients of this expansion are not constant and thus we can no longer use the 
orthonormality properties of the harmonics to extract the coefficients C[ analytically. 

Special cases can be studied numerically. For example, once the (K u , K v , K x , K y ) param- 
eters are assigned in region IV, the solution is continuously propagated into regions II and III 
without any arbitrariness, reaching the boundaries of region I. 

At each crossing point a single mode solution of region I can be numerically constructed 
and made continuous with the cited single wave solutions of region II and III. Even if they 
could be superimposed, these single mode solutions can account for a qualitative description 
of the behaviour of the Dirac field as scattered by the gravitational waves. 

In the next section the neutrino current will be analyzed separately in the four regions of 
the spacetime following this procedure. 

4 Neutrino current 

In order to study the properties of the solution and its matching across the boundary of the 
various regions of the spacetime, we analyze the associated neutrino current. 

From the Newman-Penrose formalism, the expression for the neutrino current null vector 
is the following: 

J = 2[l\F 1 \ 2 + n \F 2 \ 2 + m (F 1 F 2 *) + fh {F{ F 2 )\ . (53) 

In terms of the formalism developed so far the neutrino current density can be expressed in 
the form 

j = V=gJ = j u d u +fd v +fd x +jyd y , 

= 2{\h 2 \ 2 d u + \h 1 \ 2 d v + 2Z$le(h 1 h* 2 )d x + 2<S>5m(h 1 h* 2 )d y } , (54) 

that can easily be specialized to the various regions by identifying the correct dependence 
of hi, h 2 , 3, $. The current conservation law in the wave interaction region results in the 
simple condition 

d u j u + d v f = 0, (55) 

which trivially becomes d u j u = in region II and d v j v = in region III. Introducing the 
(closed) "projected current 1— form" 

/ = fdv - j v du , dl = (56) 

allows the use of Stokes' theorem over a bounded portion £1 of spacetime inside region I: 
J d ^ / = . In order to get information on the interaction of the neutrino current with the 
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colliding gravitational waves we examine the behaviour of the neutrino current vector limiting 
ourselves to the (u, v) plane, drawing the corresponding field lines 



du r_ f_ = \hrf ( v 

dv J" f |/ll| 2 ' 1 ' 

Along such lines the 1-form / vanishes identically by definition. 

This equation can be exactly integrated in the single wave regions, giving 

K 2 K 2 - 

% t&nv + A{v) + d (58) 



K 2 Kl 



in region II and 



K 2 K 2 ~ 

U = -JL tan W + -§ A{u) + C 2 (59) 

in region III, where C\ and C 2 are integration constants and 

A(x) = A{x) + 12. (60) 

Moreover, one can show by explicit calculation that the following remarkable property holds 
(in the single wave regions only) : the (null) neutrino current rescaled by the wavefront volume 
element 

Jr a = Va^mJ a (61) 

is geodesic, i.e. satisfies a speciality property as discussed by Wainwrightjl3j. 

The lines of the neutrino current are plotted in Figs. 2-4 for selected values of the (global) 
parameters K x and K y for a particle entering the single li-wave region from the flat spacetime 
regime, with different values of energy (i.e. for different values of K v ). In regions II and III 
the integral curves of the current density are null geodesies. 

In all the plots, the corresponding null geodesies are superimposed on the current lines. 
These are the same in the single wave region but evolve differently in the interaction region. 
Several different mode current lines are plotted, entering the collision wave region, i.e. a sort 
of mode-by-mode analytical extension of the single wave region current is attempted, since 
the general solution for hi(u, v) and h 2 {u, v) is not available here. These current lines are also 
superimposed on the geodesies. In any case, this is enough to study the general behaviour of 
the process: in the interaction region (considering the single modes for various I as elementary 
solutions of the Dirac equation) the single mode neutrino current is no longer geodesic, i.e. 
the helicity of the Dirac particle couples to the spacetime curvature; for generic values of K x 
and K y , this current shows a typical oscillatory behaviour that is more and more damped as 
one approaches the horizon; for K x = the behaviour is rather different in the sense that the 
oscillations are not as damped as in all the other cases, including the companion one K y = 0. 

It is worth noting that once again the d x direction has special properties which we have 
studied in terms of Papapetrou fields and Killing directions in a previous paper [SJ. 



5 Concluding remarks 

The Dirac equation in the background spacetime of two colliding gravitational waves with 
parallel polarization is studied in detail in the Newman-Penrose formalism. Explicit results 
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have been considered for the so called Ferrari-Ibahez degenerate solution. In this case the 
equation is separated in all the four spacetime regions, but the integration is complete in the 
single wave and flat spacetime regions only. In the collision region instead, only a partial 
analytic form of the solution can be given, involving the spin-weighted spherical harmonics. 
Continuously matching the solution across the regions' boundaries is briefly discussed too. 
The characteristics of the associated neutrino current are analysed, showing that it is geodesic 
in the single wave regions, and the results are graphically summarized (see figs. 2-6). 

Appendix. Spin-weighted spherical harmonics 

1. Analytic, implementable definition 




2. Compatibility with spherical harmonics 



Yr(e, <t>) = YT(e, <t>) 



(63) 



3. Conjugation relation 



Yr*(9, 4>) = {-l) m+s - s Y l - m {e, c/>) 



(64) 



4. Orthonormality relation 




(65) 



5. Completeness relation 




(66) 



l,m 



6. Parity relation 



(67) 



7. Clebsch-Gordan relation 




l,m,s 




(68) 
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8. Spin raising and lowering operators 
s d+ 8 Y t m = -shr< 

= VT^l Vi + s + i UiYD , 



d e (sin- s 9 s Yn + — Q d 4> (sin- s 9 s Yf 
sin 9 



sd~ S Y{ 



sm 



fl (sin 8 0,1? 



sin i 



(sin 8 9 S Y{ 



= vr+~s vi-s + i ( s -iYn . 

Harmonics with s = 1/2 and different I and m are listed in the following table. 



(69) 



(I, m) 


1/2*1 


(1/2, -1/2) 


j—e-^l 2 cos(fl/2) 


(1/2, 1/2) 


j—e^l 2 sin(0/2) 


(3/2, -3/2) 


— ^1 ie" 3 ^/ 2 cos 2 (0/2) sin(0/2) 


(3/2, -1/2) 


^y e -i<t>/2 cos (0/ 2 )(3 C os0 1) 


(3/2, 1/2) 


2 ^e^/ 2 sin(0/2) (3 cosfl + 1) 


(3/2, 3/2) 


_yi ie 3i0/2 cos ^/ 2 ) sin 2 (9/2) 
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u 




Figure 2: Several lines of the current field are plotted here, corresponding to the case K x = 1, 
K y = 1/2, 0.2 < \K V \ < 2.2, I = 1/2, where the particle is conventionally set to enter the 
single- wave zone at the point (— ir/2, 0). For v < the lines follow geodesic paths in this 
(u, v) slice, while in region I they decouple (current lines are dashed, geodesic lines are solid). 
The current clearly oscillates around u — v = const, axes, damping its oscillations as the 
horizon (the oblique line at the right) is approached. 



u 




Figure 3: Figure 2 repeated for I = 3/2, showing the influence of this parameter: it just slightly 
changes the amplitude of the oscillations, while having almost no effect on the horizon-crossing 
value of the field lines. 
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Figure 4: Figure 3 repeated with instead K y = 3/2. Neither the current lines nor the geodesies 
change qualitatively. 
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Figure 5: Figure 2 repeated for K x = 0, K y = 1/2, 0.1 < \K V \ < 1.1, I = 1/2. This choice of 
K x changes the behaviour of the current field in region I: the lines of force do not oscillate 
around directions perpendicular to the horizon (i.e. around fixed values of z) anymore and 
they are strongly different from geodesic trajectories. 
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Figure 6: This figure shows how a multiplet of current field lines forms in the transition 
between region II and region I. Here K x = 1, K y = 1/2, \K V \ = 1 and / goes from 1/2 to 9/2: 
the oscillations are damped for higher values of I. 
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-0.2 0.2 0.4 0.6 

Figure 7: Here \h\\ and \h,2\ (the widerly and finerly dotted lines, respectively) are plotted, 
with K x = 1, K y = 1/2, \K V \ = 1/8 and I = 1/2, and compared to the current (solid line). 
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